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We study, within Ginzburg-Landau theory, the responses of three-flavor superfluid quark-gluon 
plasmas to external magnetic fields and rotation, in both the color-flavor locked and isoscalar color- 
antitriplet diquark phases near the critical temperature. Fields are incorporated in the gradient 
energy arising from long wavelength distortions of the condensate, via covariant derivatives to satisfy 
local gauge symmetries associated with color and electric charge. Magnetic vortex formation, in 
response to external magnetic fields, is possible only in the isoscalar phase; in the color-flavor locked 
phase, external magnetic fields are incompletely screened by the Meissner effect. On the other hand, 
rotation of the superfluid produces vortices in the color-flavor locked phase; in the isoscalar phase, 
it produces a London gluon-photon mixed field. We estimate the coherence and Meissner lengths 
and critical magnetic fields for the two phases. 
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PACS numbers: 12.38.Mh, 26.60.+C, 97.60. Jd 



I. INTRODUCTION 



> 

Dense quark matter at low temperatures is expected to be in a BCS-paired superfluid state induced by attractive 
coupling in the color-3 channel, with the pairs having zero total angular momentum, J Because of the rich 

^-j. ■ internal degrees of freedom of quarks, two J = pairing states with condensates antisymmetric in color and flavor 
have emerged as energetically favorable candidates, a two-flavor color-antitriplet "2SC" (or particularly for u and d 

£N| . quarks, isoscalar) state, and for massless u, d, and s quarks, a color-flavor locked state The latter is the most 

stable for three flavors of massless quarks in the weak coupling limit, both at zero temperature, T, and near the 

~*^h critical temperature T c [I] . 

In this paper we study the response of BCS-paired states of quark matter to external magnetic fields and rotation, 
as would be experienced by superfluid quark matter if present in neutron stars. We find that the isoscalar paired state, 
if a Type II superconductor, behaves analogously to that of protons in a neutron star, forming vortices in response 
to the magnetic field, and forming a weak London magnetic field dominated by the gluonic component, in response 
to rotations. On the other hand, the color-flavor locked phase forms vortices in response to rotation, analogously 
to superfluid neutrons in a neutron star or to rotating superfluid 3 He, but simply tries to expel, via the Meissner 

^ . effect, a fraction of the external magnetic field p; gradients of the order parameter do not induce supercurrents, and 
magnetized vortices are topologically unstable, in contrast to the situation in the isoscalar phase. Table [| summarizes 
our main conclusions on the responses and corresponding energies. In neutron stars, the presence of the isoscalar 
phase would affect the magnetic fields of the star via the combined effects of Meissner screening, London magnetic 
fields, and possible magnetic vortices [^0|-^2[, while the color-flavor locked phase could play a role |l3| in the pinning 
and depinning of rotational vortices that give rise to pulsar glitches JlJ] . 

To investigate the nature of vortices and supercurrents induced by external magnetic fields and rotation in pairing 
states near the critical temperature, we extend the general Ginzburg-Landau theory constructed in I by adding to the 
homogeneous free energy the energy arising from long wavelength gradients of the order parameters. The gradient 
energy, which is applicable over spatial scales larger than the zero-temperature coherence length ~ T~ , is proportional 
to the superfluid baryon density n s , which we obtain explicitly in weak coupling by calculating the normal baryon 
density n n . We then incorporate external magnetic fields and rotation via covariant derivatives and transformations 
to the rotating frame, respectively. The Ginzburg-Landau free energy including color and electromagnetic gauge fields 
was analysed by Bailin and Love || and by Blaschke and Sedrakian for the isoscalar channel in the weak coupling 
limit. The present study encompasses both color-flavor locked and the isoscalar condensation of arbitrarily coupled 
systems under external magnetic fields and rotation. 

We derive the effects of uniform external magnetic fields on the condensates from the static Maxwell's equations 
in the medium, constructed by extremization of the Ginzburg-Landau free energy with respect to the gluonic and 
photonic gauge fields; the supercurrent sources include terms induced by the gradients of the order parameters, as in 
ordinary superconductors Jig]. The response of the system to uniform magnetic fields is characterized by Meissner 
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screening and free propagation of gluon- photon mixed fields, together with possible formation of magnetized vortices 
||l0|-|l2||. In clarifying the magnetic structures of the superconductor, we calculate the Ginzburg-Landau coherence 
length and the penetration depth, as well as the various critical fields. We treat all the gauge fields as averaged 
quantities and ignore fluctuations around their mean values, which cause the normal-superfluid transition to become 
first order, as clarified in Refs. |16|lT 



Magnetic fields and rotation have negligible effect on the equilibrium phase diagram near T c . Neither of these 
disturbances affects the degeneracy of the states belonging to a given phase because they do not couple with the 
orientations of the order parameters in color and flavor space. This is a contrast to the case of superfluid 3 He in 
which external magnetic fields and currents act to fix the orientations of the order parameters in spin and orbital 
space, respectively, and hence remove such degeneracy |fL8[ . In addition, the color neutrality conditions, discussed in 
I, remain unchanged, since neither magnetic fields or rotation alter color charge densities via gluon self-couplings; we 
do not take these conditions into account in this paper. While both rotation and magnetic fields yield small positive 
energy corrections that are more favorable to the isoscalar phase than the color-flavor locked phase, such corrections 
are considerably smaller than the condensation energies except in a negligibly small temperature range near T c . 

In Sec. [nj we construct the gradient energy and the supercurrents for color and flavor antisymmetric pairing with 
J = 0, and determine the superfluid baryon density and superfluid momentum density in terms of the order parameter. 



In Sec. 



Ill] we examine the responses of a color-flavor locked condensate to external magnetic fields and rotation, and in 



Sec. IV, the corresponding responses of isoscalar condensates. Our conclusions are presented in Sec. The appendix 
briefly summarizes crucial results from I. We work in units h = c = kg = 1. 

TABLE I. Responses of condensates to magnetic fields and rotation [corresponding energies] 

Magnetic fields Rotation 
Color-flavor locking Partial screening, no vortices U(l) vortices 

[Eq. (pC|)] [Eq. (^)] 

Isoscalar Parti al sc reening, possible vortices Lond on m agnetic field 



[Eq. (102)] [Eq. (125) 
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II. GRADIENT ENERGY 



To describe the effects of external magnetic fields and rotation in a quark superfiuid, we first study the gradient 
term in the Ginzburg-Landau free energy near T c . We consider a plasma composed of quarks of three massless flavors 
(uds) and three colors (RGB) at temperature T and baryon chemical potential fi, as in I, and work in terms of the 
free energy density Q = E — TS — T,i a fi ia ni a , where S is the entropy density, and the \Xi a and n^ a are the chemical 
potentials and densities of quarks of flavor i and color a. In the absence of real electromagnetic and gluonic fields the 
free density difference Ail' -' = fly' — between the homogeneous superfiuid and normal phases near T c is, up to 
fourth order in the pairing gap (see appendix), 

AQ(°) = a+Tr(0^+) F +/3+[Tr(^0 + ) F ] 2 +/3+Tr[(^0+) 2 ] F . (1) 



Here <j) + is the gap matrix for even parity, same chirality, J = quark-quark pairing (see Eq. (A.2)). The subscript 
F denotes the on-shell pairing gap evaluated at |k| = hp, the Fermi momentum of the normal state; hereafter we do 
not write the F explicitly, but it should be understood. We ignore the terms Qcn associated with the constraint of 
overall color neutrality, discussed in I, since they play no role in the ensuing discussion. 

In inhomogeneous states </>+ depends not only on the relative pair momentum, k, but also on the center-of-mass 
coordinate, r, of the pair. For wavelengths large compared with the coherence length, ~ T~ , the gradient energy 
assumes the general form to second order in spatial derivatives of the pairing gap, 

n g = ^KpTridrt+d^i) , (2) 

where di = d/dri (/ — 1,2,3). This structure follows from the invariance of the grand-canonical Hamiltonian under 
special unitary color and flavor rotation and U(l) gauge transformations of the field operators, ip — > e lv U c Uf?jj, and 
thus (4>±) Mj -» e- 2 ^{<t> ± ) cd i m {Ul) ca {Ul) db {U}) u {U}) mr 

The superfiuid mass density is related to the stiffness parameter Kt by 

Ps = ^KtTt^I) . (3) 

This relation is basically that obtained by Josephson |19| , p s = A±(rn/h) 2 \^\ 2 , for superfiuid He-II. To derive it we 
consider the situation, as in Ref. f20(| , in which the superfiuid moves uniformly with small constant velocity v s , with 
the normal fluid remaining at rest. The momentum per baryon in the flow is yuv s . In such a situation the phase factor 
of the gap in the lab frame transforms by <p+ — ► e -lP ' r </>+, where P is the total pair momentum relative to its value 
in the lab frame. Since each pair carries a baryon number 2/3, 

P = \^s- (4) 

Under the transformation cf> + — > e~ lP ' r (f> + the free energy density transforms by 

n - n - l -K T P ■ Tr(0 + V^ + - ^V0+) + ^ T P 2 Tr(^+^) . (5) 

Comparing this result with the expected transformation of the free energy density, 

n -»• O + g s • v s + ^p s v 2 s , (6) 

and using Eq. (Q), we derive Eq. (^), and in addition relate the momentum density of the condensate, g s , to the order 
parameter by 

K. = - l -K T ^{4>+V4>l - 4>l^4>+) - Ms , (7) 

where j s is the superfiuid baryon current. [The minus sign originates from the conventional definition of the gap here 
in terms of the adjoint spinors rather than the spinors, as in condensed matter physics.] 
In the limit of weak coupling near T c , the stiffness parameter is given by 

7C(3)n h 

Kt ~ WMfr ' (8) 
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where n& is the baryon density, and £(3) = 1.202 . . .. This result has been derived by Giannakis and Ren plj using 
finite temperature diagrammatic perturbation theory in the normal phase for single gluon exchange interaction, and 
by Bailin and Love Q for a short-range pairing interaction. Here we present a simple physical derivation of Kt from 
Eq. (H), by evaluating the supcrfluid mass density p s near T c in weak coupling. We consider the system, as in Leggett's 
argument for superfluid 3 Hc ]l8| , to be contained in an infinitely long tube which moves uniformly with velocity v 
(|v| <C c) along the symmetry axis. We assume that the normal component is in equilibrium with the moving wall 
and that the superfluid component stays at rest. Then the baryon current density j„ of the normal component is 
related to the normal baryon density, n n = nt, — n s , by 

jn = n n v , (9) 

and is given in weak coupling in terms of the quasiparticle distribution function, /k, by 

d 3 k k „ 

-Tr/ k . (10) 



(2tt) 3 (i 

In the present situation, /k is the equilibrium distribution evaluated in the frame of the moving walls: 

/ k = f(e(k) - k ■ v) = e[e(k) _ k 1 v]/r + 1 ; (11) 

here 

£ (k) = [(|k|- M /3) 2 + ^(k)0+(k)] 1 /2 , (12) 
with 0+(k) = 0+(e(k), k), is the energy matrix for quark quasiparticles [I]. Expanding to linear order in v we find 

d 3 k k 2 



n n = - 



J (2^^T Tr[/(£(k))(1 ~ /(£(k)))] ' (13) 



The integration is dominated by momenta k close to the Fermi surface, where the gap can be taken to be constant. 
Expanding near T c to second order in the gap we thus obtain 



1 - IS? T *^> 



n b , (14) 



so that 



n s = p s /fi = -I^L T r(<^0 + )n & , (15) 

from which Eq. (|) follows. 

Real color and electromagnetic gauge fields enter the gradient energy (^) via covariant derivatives. Let us first 
consider color SU (3) local gauge transformations of the quark spinors via 

U c (r) =exp[zA> Q (r)/2] , (16) 

where the ip a are the local phase angles and the A a the Gell-Mann matrices. The gap matrix <fi+(k; r) transforms as 
(4>+)abij — > (Ul) ca (U^)db(4'+)cdij, and the corresponding covariant derivative is 

D t 4>+ = M+ + Ug[{\ a Y^+ + 0+A"]A« , (17) 

where g is the qcd coupling constant, and the A" are the color gauge fields, which transform as Af — > Af + 
g~ 1 d[Lp a + fa/3-yAfipj. With extension of the covariant derivative ( p~7[ ) to include the usual U(l) electromagnetic 
gauge transformations, which multiply the quark spinors ip a i by 

Ui(r) = exp[igi¥? e (r)] , (18) 

with local phase angle <p e , the full covariant derivative is then pll 

Di4> + = + ^g[(X a )*^ + + <f) + \ a ]Af + leQ^+Ai , (19) 
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where A is the electromagnetic gauge field, e is unit of the electric charge, and Q a bij — S a b(qi + qj) is the electric 
charge matrix of the pair. The full gradient energy is then 



n g = -^ T Tr[A0+(A0+) t ] 



(20) 



with D\ given by Eq. (|l^). This form is invariant under transformations (4>+) a bij — * U*U* (U£) ca (U^)db(4>+)cdij ■ 

We now write the free energy densities, f2„ and f2 s , of the normal and superfluid phases near T c in the presence of 
an external magnetic field, H cxt , neglecting the response of the normal component to H cxt . First, 



a 



|H„ 



(21) 



where F n is the free energy of the normal state at H cxt = 0. To construct the free energy density of the superfluid 
state, we add to the homogeneous free energy (|l|), the gradient energy density (f20|), and the energy densities of the 
color and electromagnetic fields induced by H oxt , to obtain 



n s = F n + a +Tr(<^+) + /3+[Tr(^ + )] 2 + /3 2 + Tr[(^0 + ) 2 

1 



-^ T Tr[A<MA0+) t ] +\G? m Gf m 



:Fi m Fi 
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where 



and 



Go 
h. 



-dtA^ + d^-gf^AfAl 



lm 



-diAr, 



d m Ai 



are the field tensors. 

Extremization of J d 3 rA£l, where AJ7 = f2 s — Q n , with respect to (fy, yields the gap equation 
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K T D l (D^ + )a+cf, + + 2/3+[Tr(0U + )]«^ + + 2/3+0+^0+ = . 



(22) 

(23) 
(24) 

(25) 



The field equations for the macroscopic fields A a and A, found by extremizing J G? 3 rASl with respect to the fields, 
are: 



9fai3'yA^ n G^ nl 

l -K T glm {Tr [((A Q )*^+ + 0+A") 1 %</>+] } 
K T g 2 A?Re {Tr [((\ a )*<t>+ + 0+A") ((^ /3 )*0+ + ^+A /3 ) t ] } 
KrgeAiRe {Tr [((A q )>+ + 0+A") Q^ + ] } = J? 



(26) 



and 



d m F ml = -K T elm[Tr{Q4>ld^ + )] - K T e 2 ' A{^{Q^ 
1 

~2 J 



^rgeAfRe {Tr ((A«)*0+ + + A«) f ] } = Ji , 



(27) 



where J Q and J are the color and electric supercurrent densities, respectively. The Maxwell equations ([26|) and (|27| ) 
with uniform external fields determine the structures of the supercurrents and vortices, as we discuss in Sees. Ill and 



IV 



For later analysis it is useful to write down the solution of these equations in the presence of an externally applied 
electromagnetic transverse field, A oxt (r) = e lq r A cxt (q), where q • A cxt (q) = 0. The supercurrents, J and J Q , induced 
by this potential in turn induce fields, A; n d and A^ ld . We assume that the external potential A ext is sufficiently weak 
and nearly uniform (q — > 0) that the order parameter remains essentially homogeneous and the self-couplings of the 
induced gluon fields are negligible. Up to linear order in A cxt , the induced fields may be calculated from Eqs. ( |2^ ) 
and (p?|), with Gf m and Fi m including the induced fields. The supercurrents J and J" depend on the total fields, 
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\a _ a a 

— ^ind ' 



(28) 



which are given by 



A(q) 



x -Re 
4 



xRe 



Tr(Q<M(AT0+ + 0+A a )' 
Tr (((A /3 )>+ + </!>+A /3 ) Q^)] A cxt (q) 



and 



A°(q) = -^[(e i ) ]a/3— — 



Re 



Tr 



(((A")V+ + Q^)] A(q) , 



(29) 



(30) 



where the transverse color dielectric function, 



_T r 

£ q/3 — °a/3 



+ ^Re {Tr [((A Q )>+ + + A Q ) ((A^)*0 + + 0+A^) f ] } , 



(31) 



describes the response of the current of color a to the field of color f3. In deriving this equation from Eq. (34) of II, 
we have written the long wavelength limit of the transverse screened correlation function j^T (k, 0) in terms of the 
gap as 



X?(k - 0,0) = \K T g 2 Re {Tr [((A«)*0+ + 0+A") ((A^)>+ + 0+A^ 1 ] } 



(32) 



In the succeeding sections we analyze the structures of inhomogeneous color-flavor locked and isoscalar condensates. 
These condensates are antisymmetric in color and flavor, with a gap 



4>+ — £ijh e abc{d c )h : 



(33) 



where (d c )h is a complex field representing the pairing gap between two quarks in state [(\ab) — \ba))/y2] x — 
\ji))/\/2] with a ^ b ^ c and i ^ j ^ h. With (|33]), Eqs. (22) and (|2l]) for fl s and Q„ lead to a thermodynamic 
potential density difference near T c 



Att = a\ + (01 + /3 2 T)A 2 

+K t Tt[DMD 1 ^} + -G? m G? m + -F lm F lm - i|H ext | 



(34) 



where <S> a bi = £abc(d c )i. The homogeneous part is that in Eq. ( A8). 

We turn now to the properties of the color-flavor locked and isoscalar condensates, with respective gap matrices 
( [Ail ) and (|Al|), in the presence of external magnetic fields and rotation. 



III. COLOR-FLAVOR LOCKING 



In the color-flavor locked phase, vortices associated with the response of the system to magnetic fields are topolog- 
ically unstable. Rather, the system acts as an imperfect diamagnet. On the other hand, rotation of the color-flavor 
locked phase produces vortices associated with the baryon U(l) symmetry breaking, analogous to the vortices in 
rotating chargeless superfluids such as liquid He II Q, superfluid 3 He [Q, Bose-Einstein condensates of alkali atoms 
p5|, and neutron superfluids P,E6|. 



We consider the color-flavor locked phase described by the gap (A14) 



(<t>+)abij = K A {5ai5bj ~ S a jS bi ) . 

In the absence of magnetic and color magnetic fields the gap equation becomes 



(35) 



G 



an A + 6P cfl \k a \ 2 k a - 2K t V 2 k a = , (36) 

where /?cfl = Pi +02/3- In the London limit, in which the scale L of the spatial variation is much larger than the 
Ginzburg-Landau coherence length, 

?cfl - (2K T /\a\) 1/2 , (37) 

the relative variation of the magnitude of n A is of order Ccfl/^ 2 smaller than the relative variation of the phase of 
k a and can be neglected. To see this we write n A = e lipA \n A \. In the homogeneous limit, ip A can be taken to be zero, 
and the magnitude of the gap is given by 

a + G^cFLM 2 = . (38) 

Thus in an inhomogeneous situation the variation of the magnitude of the gap is related to the variation of its phase 
by 

T^T = S&x^VV^ . (39) 

\K A \ 2 

Equation (j38|) implies the homogeneous phase condensation energy 

a 2 

-AQ= -2 — . (40) 



A. Response to magnetic fields 

To determine the response of the condensate to an external electromagnetic vector potential we substitute the gap 
( |35| ) into the field equations ( p6| ) and (]27]), and neglecting gradients of the magnitude of k a , we find the color and 
electromagnetic supercurrent densities: 

J" = -2K T g 2 \ KA \ 2 A a , a ^3, 8, (41) 
J 3 = -2K T g\ KA \ 2 (gA 3 + eA) , (42) 

-2K T g\ KA \ 2 ( ff A 8 + ^=A) , (4:-!) 



-2K T e\ KA \ 



V3 



-eA + q [A 3 + ^A 8 

3 y \ V3 



(44) 



Note that the current terms induced by the gradients of the phase of k a vanish due to the isotropy of the condensate 
in color and flavor space, a feature implicit in Refs. p] , p7| . This absence of currents induced by the gradients of the 
phase of k a underlies the response of the color-flavor locked condensate to magnetic fields and rotation, as we consider 
in the present and following subsections. External disturbances such as magnetic fields and rotation do not c hang e 
such isotropy and hence do not break the degeneracy in energy of the states with order parameters of the form ( |Al3j ). 

The currents associated with color charge a — 3,8 and electric charge respond to the electromagnetic gauge field; 
in turn the electromagnetic current responds to the gauge fields A 3 and A 8 . As in Ref. it is convenient to take 
linear combinations of the supercurrents and the gauge fields to diagonalize the response. The combinations 

_ V3.gJ - V3eJ 3 - eJ 8 ^ VSgA - V3eA 3 - eA 8 
= 3y/2g 3 ' = 3\/2 ff 3 ' 

3 ^ 3, 9 J 3 + 4eJ + V3. 9 J 8 A3 _ 3gA 3 + 4eA + ^gA 8 
6V2 53 ' &V2g 3 
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r , s _V3J 8 -J 3 A8 _V3A 8 -A 3 



J=^— 2 > 2 (4?) 

satisfy the diagonalized equations 

J = , (48) 



J 3 = -12K T g 2 \n A \ 2 A 3 , (49) 

and 

J 8 = -2K T g 2 \n A \ 2 A\ (50) 

where 

.93 = ^V / 3ff 2 + 4e 2 (51) 

is the coupling constant associated with the field A. 3 . Since Eqs. (pd|), (|49|), and ( |50| ) imply that the corresponding 
fields have a magnetic mass, we see the property originally shown by Alford, Rajagopal, and Wilczek Q, that eight 
of the nine gauge fields are Meissner screened, while the field A remains free. 

The response of a homogeneous color-flavor locked condensate near T c to an electromagnetic field is twofold. 
Substituting the supercurrents (|4"l])-(f44|) into the linear response (|2^) and (|30"|) we obtain 



and thus 



q 2 + 2K T g 2 \n A \ 2 A 3_./o A 8 2K T ge\ KA \ 2 

A -q 2 + l2K T gl\n A \^ ' ~~ _ <? 2 + 12^ T5 2 M 2 Acxt ' (52j 



3 A A 3 = ^ 2 , * 2| |2 A ext , A 8 = 0. (53) 



V6g 3 ' 3. 93 g 2 + 12X T g 3 2 l^l 2 

The term A 3 , which is proportional to A cxt , is screened by the Meissner effect within a penetration depth 

Acfl = * ■ (54) 

2y/3K T 93\K A \ 

[However, see note (8).] The mixed field, A, defined by Eq. (45), also proportional to A ex t, propagates freely in the 
superconductor. The behavior of the field A 3 is comparable to that in a Type I superconductor, independent of 
the relative magnitudes of the coherence length and the penetration depth. Due to the absence of phase-gradient 
induced current terms, an external field does not lead to vortex lines; a vortex placed inside the condensate would be 
topologically unstable. We note that since the structure factor / a 38 vanishes, the gluon self-coupling term included 
in the field strength tensor (|23|) does not appear, allowing us to generalize expression (|5^ ) straightforwardly to the 
case of finite external fields. 

Under the external field H cx t = V x A oxt , the free energy difference near T c , Eq. (|34|), reduces to 

A SI = 3a\ KA \ 2 + 9/3 C flM 4 + 6K t (\Wk a \ 2 + 9 2 \A 3 \ 2 \k a \ 2 ) + X - (\B 3 \ 2 + \B\ 2 - |H cxt | 2 ) , (55) 

where 

B = V x A, B 3 eee V x A 3 (56) 



are the Abelian field strengths. The gap equation (cf. (|36|)) follows directly from Eq. (p5[). 

The system remains a superfluid for magnetic fields H cxt below a critical field H c . For H cxt < H c , the equilibrium 
system expels the part of the external field that leads to B 3 . An external field greater than H c entirely penetrates the 
system, and the system becomes normal. The favored phase is the one that minimizes the Gibbs free energy density, 
G = SI — H cxt • B, at constant external field, H cxt . Since B 3 is screened out in the color-flavor locked phase in bulk, 



a~7, 1" 77 1^| 2 ~ H cx t ■ B — F n — —= ~Sr? 

4/?cfl 2' 4/3 C FL 12fl§ 



G s = F n - ^— + -\B\ 2 - H oxt -B = F n - ^IH^I 2 . (57) 
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The normal phase Gibbs free energy density is 



G n — F n — — | H C xt | 2 , (58) 



so that the difference of Gibbs free energy densities of the superfLuid and normal phases is 



with the magnetic part 



AG = ——= + AG mag (59) 

4pcFL 



AG mag — - |2?| 2 — H cxt ■ B + -|H 0Xt | 2 



c 



2 



2|H cxt | 2 . (60) 



9ff 3 

This free energy difference vanishes at the critical field 

2eV/3 C FL ^ e sCFL AcFL 

The coh erence length (|37|), the penetration depth and the critical field ( |6l| ) can be estimated from Eqs. 



flAq)-( A6) and (p[) in the limit of weak coupling as 



, 1 \ /300 McV\ / T . 



-1/2 



ff c * 1.8 x 10 19 ff3 ( —5—) f (l-£)G. (04, 



and 

,100 MeV 7 V300 MeV ) \ T ( 

In these expressions, T c and \x are normalized by their typical values at lower densities. Although strictly valid near 
T c , these estimates of the characteristic lengths and critical field can be extrapolated to T -C T c with the given 
temperature dependence in terms of 1 — T/T c without modification, an approximation which should be accurate to 
within a factor of two, as in the case of ordinary superconductors pl|. The estimate ( |64] ) indicates that a color-flavor 
locked condensate, if present in a neutron star core, would not be destroyed by ambient neutron star magnetic fields. 



B. Response to rotation 



We turn to analyzing the response of a color-flavor locked homogeneous condensate near T c to rotation. For 
simplicity, we consider the condensate to be in thermodynamic equilibrium in a cylindrical vessel rotating with 
constant angular velocity a;. In this situation, the normal component corotates with the vessel as a rigid body. The 
superfluid component, however, behaves independently; to examine its properties, we transform to the rotating frame 
in which the container walls are at rest, assuming, as in a neutron star, that the velocities associated with rotation 
to be much smaller than the speed of light and describing the flow properties nonrelativistically. 

The free energy density in the rotating frame is given by £1 — u) ■ L, where L is the angular momentum density, 
given in terms of the order parameter by 

L = r x g s , (65) 

with gg given by Eq. (0). We neglect here the normal component contribution to the angular momentum. Using 
Eq. (|35|), we find that the superfluid mass density, (^), is 







Ps = ^K T p 2 \ KA \ 2 , (66) 

and 

g s = UK t ^{k* a Vk a - kaV4) • (67) 

Thus in the rotating frame, the gradient energy included in Eq. ( j55[) becomes 

n g = 6K t {\Vk a \ 2 ~ (2ili/3)u> • [r x (k\Vk a - k a Vk* a )] + g£\A 3 \ 2 \K A \ 2 } , (68) 

with A. 3 given by Eq. (|46|). Note that the gradient term is independent of the vector potentials; thus rotation produces 
a lattice of vortex lines, with the coarse grained s uperfl ow pattern simulating corotation of the condensate. In contrast, 
in an isoscalar condensate, as we discuss in Sec. IVB, a rotating condensate produces a London magnetic field rather 
than rotational vortices. We assume in the following zero external magnetic field. 
A singly quantized rotational vortex centered on the cylinder axis has the structure, 

k a = e-*+\K A \ , (69) 

where <j> is the azimuthal angle around the vortex line. The corresponding baryon current density, defined by Eq. (Q), 
is 

j s = 8K t \k a \ 2 V4> = ^ V^ = n s v s , (70) 
2 ll 

where we relate k a to n s via Eq. (|66|). Integrating v s around a closed loop surrounding the vortex line, we have 

dt ■ v s = 2tt^- ; (71) 
2\i 

thus an individual singly quantized vortex has circulation 37r//i. Specializing to a circular contour of radius r centered 
on the line, we see that the vortex velocity, in the azimuthal direction, is given by 

|v s | = ^-, forr>£ C FL. (72) 

In addition the energy per unit length of an isolated vortex is 

9irn s ( R \ 

T L = —^\n , (73) 



4//. 



,CFL 



where R is the container radius. The onset of vorticity occurs at the critical rotation rate u> c i, where the energy in 
the rotating frame Tl/ttR 2 — ui ■ L first vanishes; thus 

At finite rotation speeds 3> uj c i the system forms a triangular lattice of singly quantized vortices. The overall 
flow of the condensate is essentially that of a solid body at angular velocity u, with velocity, Eq. (f72j), close to the 
vortex cores. This situation is similar to the case of superfluid neutron vortices in a rotating neutron star where 
the quantization of circulation is 7r//i 0. The net circulation around the cylinder is 3irN v /[i, where N v is the total 
number of the vortices. Since at the boundary, |v s | ~ coR, the circulation is 2irR 2 tu. Thus 

2 , 

V Prot J V300 MeV J V 10 km J ' V ' 

where in the latter equation the quantities P rot (= 2tt/o;), /x/3, and R are normalized to characteristic neutron star 
scales. For these values, the intervortex spacing, ~ (nR 2 /Ny) 1 / 2 , is much larger than the Ginzburg-Landau coherence 
length £cfl, except immediately near T c . 

The total angular momentum, L to t, of a vortex line, the integral of (j,t X j s , is 
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Ltot = -N s z 



(76) 



where z is the direction of the vortex axis, and N s is the total superfluid baryon number. As in a neutron superfluid 
in a neutron star, a lattice of vortices in a color-flavor locked superfluid is an effective store of angular momentum. 

To calculate the superfiow energy density, O rot , of a vortex lattice in the rotating frame we note that in the 
logarithmic integral for the kinetic energy per line, the radial integration is cutoff at the vortex lattice constant a, so 
that the line energy is given by Eq. (73) only with R in the logarithm replaced by a = (2y/3ir/ fiui) 1 /' 2 ; thus in the 
rotating frame, 



i 



-p s (w x r) 2 , 



(77) 



where the latter term is the free energy in the rotating frame of the (coarse grained) corotation of the fluid, and 

3n 



' 111 



CFL 



(78) 



is the local vortex contribution. In weak coupling, ASl rot can be estimated from Eqs. (A5)-(A6) and (|8j) as 

3 



Aft rot ~3x 1(T 18 



( M/3 



V300 MeVy \P T 



1 ms 



1- — I In 



CcFL 



MeV fm 



-3 



(79) 



The upper critical velocity, w c2 , at which rotation destroys superfluidity, occurs where the intervortex spacing ap- 
proaches the coherence length, an impossibly high velocity to be relevant for neutron stars. 



IV. ISOSCALAR, COLOR- ANTITRIPLET CHANNEL 



The response of the isoscalar, color-antitriplet phase to magnetic fields and rotation is quite different from that of the 
color-flavor locked phase. Magnetic fields lead to topologically stable vortices, as in laboratory Type II superconduc- 
tors, while rotation of the condensate induces a London magnetic field p8| , p9fl , rather than rotational vortices. Since 
these external disturbances do not distinguish any direction in color space, the color structure of the order parameters 
is unaffected. We specialize, without loss of generality, to a gap matrix having the specific color orientation: 

(daJise^ldlWi. , (80) 

where B denotes the blue and s the strange directions; then 

{<P+)abij =ey,e aU} e , '* ,0 |d|. (81) 

We concentrate on the London limit in which the spatial gradient of the magnitude of the gap |d| can be neglected (see 
Eq. (|3^)) and that of the phase tpo varies very slowly. This approximation is sufficient to describe the supercurrent 
structure outside the vortex cores. 



A. Response to magnetic fields 



Substitution of the gap ( |8l| ) into Eqs. ( pq ) and ( p7| ) leads to the relations between the supercurrent densities and 
the vector potentials, 

J" = 0, a =1,2, 3, (82) 



J 4 = -K T g\d\ 2 (gA A - 2ImV4) , (83) 
J a = -K T g 2 \d\ 2 A a , a = 5, 7, (84) 
J 6 = -K T g\d\ 2 (gA 6 - 2ImVd G ) , (85) 
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^-^rflldr ( V<^o + ^A 8 + -A ) , (Mi) 



(87) 



V3g 



where d = d/|d|. 

As in the color-flavor locked case, it is useful to diagonalize these relations via the linear combinations ]23j 

J^* 3 -* 3 *, A=^ A - eA \ (88) 
358 3ff 8 

3#8 3g 8 

where 

J = (90) 

and 

J 8 = -4K T g 8 \d\ 2 (V^ + g s A 8 ) , (91) 

with 

.98 = Jv / 3^+^ (92) 

the coupling constant associated with the field A. 8 . Thus five of the nine gauge fields are Meissner screened, as 
expected from the results for color Meissner effects on virtual gluons mediating quark-quark interactions (see, e.g., 



Ref. 1 30 1 ) . Expressions (90) and (91) indicate that the photonic gauge field, A, modified by a gluonic component 
of charge a = 8, is a free field, whereas the field A 8 couples to the gradient of the phase <po. This latter coupling- 
underlies the response of the isoscalar homogeneous condensate to external magnetic fields and rotation. 

Equations (^9|) and ( |30| ) determine the linear response of the system to weak external magnetic fields. For the 
order parameter (|Sl|), an external electromagnetic vector potential leads to terms in these equations proportional to 
Xg B = A33, which in the standard representation is non-zero only for a = 8. Thus only a field A 8 is produced; the 
fields A" remain zero for a = 1-7. The absence of the a = 4 and 6 gauge fields implies via Eqs. (^3j) and (|5|) that in 
equilibrium external magnetic fields do not lead to non-zero gradients ImVdR and IrnVdc. For the electromagnetic 
and a — 8 gauge fields we obtain 

_ g 2 + (4/3)g T ,g 2 |d| 2 (4/3V3)K T ge\d\\ 
A q2 + 4K T gl\d\2 Acxt ' q 2 + 4K T g 2 \d\^ Acxt ' 



or, equivalcntly, 



9 * a* e I 2 



A V3 58 AcXt ' A 3. 9 8 q 2 + 4X T .g 2 |d| 2 Acxt ■ (94) 



We note two features similar to the case of the color-flavor locked condensate. First, the modified field A is subject 
to the Meissner effect with a penetration depth (!(]] 

A IS = )=— ; (95) 



the modified field A is free. Second, the gluon self-coupling term in the field strength tensor ( J23| ) is absent. This fact 
allows the extension of Eq. (|94|) to finite external fields. 

The free energy difference near T c , Eq. (|34|), may now be written in terms of A and A 8 as 

An = a\d\ 2 + /3 IS |d| 4 + 2KtM + *g s At)d\ 2 + ±\B 8 \ 2 + \\B\ 2 - i|H cxt | 2 , (96) 
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where 



B = V xA, B s = V xA s 



(97) 



are the Abelian field strengths, and (3is — Pi +Jh- As in ordinary superconductors |l5[ , we identify the Ginzburg- 
Landau coherence length, £iS; from expression rt9q) as 



6s = (2K T /\a\) 



1/2 



(98) 



a result identical with the color-flavor locked expression (|37|) . 

The thermodynamic critical field, H c , associated with B = V x A, is the field at which the Gibbs free energy of 
the normal state drops to that of the superconducting state. Since B s is screened out in the isoscalar phase in bulk, 



Gs — F n 

where we have used the identity 



IS 



\\B\ 



B = F„ 



4/? 



is 



—2 |H ex t| 2 
6<?8 



(99) 



■J—B s 



9 



derived from Eqs. (|8|), (|89|) and the first of | 
normal phase is thus 



H ext , (100) 
The difference in Gibbs free energy densities of the paired and 



AG 



a 



Aft_ 



IS 



(101) 



with the magnetic part 



AG,, 



B 



2 I H ext | 



|Ho 



Then, 



H r = 



e ^2As \/2e6sAis ' 



(102) 



(103) 



The presence of the gauge field, A , in the derivative term in the free energy indicates that the system can support 
magnetic vortices, as in Type II laboratory superconductors. Whether or not vortices appear depends, as we discuss 
below, on the relative size of the screening and coherence lengths. Let us first analyze the structure outside the core 
of a single vortex, whose core, of radius ~ £isj w ^ take to be aligned along the z-axis. Then the order parameter of 
the vortex is given by Eq. (p0|), with ipo = —(f>, where tfi is, as before, the azimuthal angle around the line. The field 
equation for S 8 , derived from extremization of the integral of Eq. (|9^) over the system with respect to A s at fixed 
A, is 



8 



1 



Af s V x (V x B s ) = —V x 



2tt 

38 



5(x)S(y)z 



(104) 



in the London limit in which the spatial gradient of the ma gnitu de of the gap |d| is negligible. The condition 
of quantization of magnetic flux, found from integrating Eq. (104) over the interior of a contour in the x-y plane 
surrounding the line, is 



£d£-{A 8 



98 



(105) 



The flux quantum, 08, is thus 2it/g%; note that in the absence of color coupling (g — 0) the flux quantum reduces to 
6-7r/e, which is 2ir divided by the net charge, 2/3 — 1/3, of the pair. 
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As in ordinary Type II superconductors JL5|], it is straightforward to solve the London equation (104). The current 



density , which flows in the azimuthal direction around the line, has a magnitude 



2ir\j s r ' 



for £i S < r < Ais , 



1/2 



T^ + 7^)e- r/Als , for r» Ais 
Ais 2r , 



(106) 



and the current is screened by the Meiss ner effect far from the vortex core (r 3> Ais); there the gauge field A. s 
alone fulfills the quantization condition (105). The vortex line energy per unit length, the sum of the magnetic and 
flow energies, is 



Tr = 



4^ Uis 



(107) 



We proceed to classify the superconducting properties according to whether or not the isoscalar condensate allows 
magnetic vortices associated with the field A. 8 to form, i.e., whether it is Type II or Type I. Blaschke and Sedrakian 
JTo| addressed this problem in the weak coupling Ginzburg-Landau formalism. We give the argument here for arbitrary 
coupling. The basic idea is to calculate the energy per unit area, <r s , needed to form a planar surface separating the 
normal and superconducting material. At the thermodynamic critical field H c , applied parallel to the surface, the 
surface is in mechanical equilibrium. For a surface perpendicular to the z-axis, the surface energy a s per unit area 
may be written as the integral over z of the difference of the total Gibbs free energy density, A£l(z)\H cxt ^H c + F n + 
\H"l — H c \B(z)\, and the value F n — \H"l for the limiting case in which the thickness of the interface vanishes (Ais — * 
and £is — > 0); we find 



/oo 
dz [An(z)\ H ^ Hc - H c (\B(z)\ - H c )} 
-oo 

dz{a|d(^)| 2 +As|d(z)| 4 



+2K T \[d l +ig s A 8 (z)}d(z) 

-\ 2 ' 



1 

+ 2 

6s 



\B 8 (z)\ - - — H, 



358 



-H, 



3gs 



(1 - V2k is ) , for Kjs ~ 1/V2 , 



(108) 



where the Ginzburg-Landau parameter «is is defined by 

«is = Ais /6s 



(109) 

with Ais an d Cis give n by Eqs. (|9^) and (|98|). The calculation of the prefactor of a s must be done numerically. The 
structure of Eq. (|l08| ) is identical to that of ordinary superconductors [[l5[ with the surface energy changing sign from 
positive to negative as k\s goes through l/v2; thus we are led to the same criterion, namely that the system is Type 
I, with no vortex formation, for a s > 0, and Type II, with an Abrikosov-Schubnikov vortex phase, for a s < 0. 

While in a Type I superconductor, superconductivity is destroyed by an external uniform magnetic field H cxt greater 
than iJ c , in a Type II superconductor, vortices form at a lower critical field H c \ < H Cl and increase in density up to 
an upper critical field H C 2 (> H c ), where the system turns normal. At the field H C 2, magnetic vortex cores essentially 
fuse. At the lower critical field H c \, the energy gain from penetration of the magnetic field just compensates for the 
energy loss due to the line energy Tl, Eq. (107). The energy gain per unit volume from field penetration, given in 
terms of the magnetic induction B 8 and the corresponding magnetic field H s — (e/3gs)H cxt (cf. the s econ d part of 
Eq. (|39j)), is E mag = —B s Ti 8 . Integrating over the volume including the vortex, we find, using Eq. (105), a total 
energy gain from field penetration, (e/3g8)i? C xt08- Equating this term to Tl we find that 



3 /Ais 
2 e A? s m Uis 



(110) 



Just below H C 2 the pairing gap has a tiny magnitude but strong spatial variation, and is determined by a linearized 
and inhomogeneous version of the gap equation (p5[), as first discussed by Abrikosov |31|. From the condition for the 
existence of nontrivial solutions to this equation, we obtain the familiar form 
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Hc2 = ■ 

efts 

At high densities, where the system is weakly coupled, kjs reduces to 



Kis = 18tH 



V7C(3) 



1/2 



Tc_ 



(in) 



(112) 



Since T c ~ /ig~ 5 e~ 37I ' 2 / v ' 2 ~ 9 ]3^j , kjs is far smaller than unity, so that the system is Type I near T c . It is not out of 
the question that at low densities where T c can be ~ 0.1 \i and g$ ~ 1, that K^g can become larger than l/v2 and 
the system can become Type II |^o| , expelling the field B 8 completely for ff ex t < ffcii while allowing the external 
mag netic field i? e xt to penetrate freely in the form of the photon-gluon mixed field 3 of strength (g / 1 ^fSg§)H cy ± [Eq. 

©]■ _ _ " ' _ 

The magnitude of the coherence length fts, fl9q), is the same as £cfl in Eq. (|62j), while the penetration depth pq ) 

can be estimated from 



A 



is 



1.5 



/ 300 MeV \ 
g a I V M/3 / 



-1/2 



fm 



To estimate the magnitude of the critical fields ( |103| ), ( |11 
H c ~ 3.6 x 10 19 V 



and (111) near T c we use ( A5 )— ( A.6 ) and 
M/3 



300 McV 



T 



G , 



(113) 

to find 

(114) 



H cl ~ 7.9 x 10 18 ( ^ 



/ /i/3 



V300 MeV 



and 



1 + ln 



V3 T c 



300 MeV 



38 100 MeV n/3 



1- — ) c; 



(115) 



c2 



2.9 x 10 2 



100 MeV 



1 



(116) 



Extrapolation of expressions (114)— (116) to low densities and temperatures indicates that the critical fields H c , H c i, 
and H C 2 are several orders of magnitude larger than canonical neutron star surface fields ~ 10 12 G. However, to assess 
the actual situations possible in neutron stars, one must take into account the history of the expulsion of the magnetic 
field and the possibility of freezing in of the magnetic field; see note pj . 



B. Response to rotation 

To determine the equilibrium properties of a rotating isoscalar condensate, we assume, as in the color-flavor locked 
discussion, that it is in a cylinder rotating at angular velocity u). We first note that Eqs. (||) and ( [sif ) imply that the 
superfluid mass density of an isoscalar condensate is 

1fi 

Ps = y W|d| 2 ; (117) 
in addition, the momentum carried by the condensate is 

G;.), = ^tmK • (iBi - gsAf)d + [(-idi - g s Af)d*} ■ d} . (118) 



Thus, in the rotating frame, the gradient energy in Eq. (96) is 

tt a = 2K T \(di + ig s Af)d\ 2 - *K Tf i(w x r),{d* • (idi - g & A^)d + [(-idi - g s Af)d*} ■ d} , (119) 
where Af is given by Eq. 
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From expression (119) we find that in equilibrium, rotation induces a London magnetic field associated with the 
gluon-photon mixed potential A 8 , rather than generating vortex lines as in a color- flavor locked condensate. This 
response is due to the presence of charge g$, and is thus similar to the case of ordinary rotating superconductors 
|)28| , [29f and proton superconductors in rotating neutron stars [^|. Rotational vortices are topologically unstable; a 
rotational vortex of the form (173) in the system can be unwound by bringing the velocity to zero simultaneously 
generating a gauge field A 8 to preserving the quantization condition ( 105 ). To derive the London field we assume 
that the system has zero external magnetic field. In a rotating superconductor the local current JT 8 is given, deeper 
in than a penetration depth from the surface, by 



^8 3 

J = ^s"^ x r ) 

since the pairs, of charge gs and baryon number 2/3, corotate with the vessel, 
quantization condition on the vorticity becomes 



(120) 



In the absence of a vortex, the 



j> dl ■ (A 8 + Xj s J s ) = 



(121) 



[cf. ( 105 ) for a contour surrounding a quantized vortex]. Taking a circular contour of fixed distance r from the rotation 
axis, we find 



,11 ■ A'' = / d 2 rBl = -P^cur 2 . 

3# 8 



(122) 



Thus B, is constant and 



6 = -tt-w, *4 = -^-w x r 
3.98 og s 



(123) 



The magnetic field B remains zero. 

Since d is a constant in space, the gradient energy (119) reduces to — \p s {u x r) 2 so that the total free energy 
density associated with rotation, measured in the rotating frame, is 



n rot = Afi ro t - -p 8 (a>xr) 



(124) 



where 



Afi rc 



(125) 



is the energy density of the London field. 

We make several observations about the London field (123). First, it reduces to the conventional London field of a 
nonrelativistic electron superconductor with the replacements g%/2 — > —e and /u/3 — > m e . Second, since B = 0, we 
have |B 8 | = (\/3.9/e)|B|, which implies that the London field is primarily composed of the color field of charge a = 8. 
The order of magnitude of the London field is 



IB 8 



0.15 




(126) 



so that at low densities relevant to neutron star interiors, the London field is thoroughly negligible compared with the 
typical magnitude (~ 10 12 G) of the neutron star surface magnetic fields. Third, in the presence of a uniform external 
field £f ox t < H c \, we have B — (5/v / 3g8)H ex t, while the London field ( |126| ) is unaffected. For if ext > H c x, where 
vortices form, we find that the total field B remains the same, while the total field B 8 is the sum of that associated 
with the vortices, plus the London field, as long as the vortex separation is large compared with the penetration 
depth, Ais- Fourth, the London field increases with increasing density, because of the power-law increase in /i/3 and 
the logarithmic decrease in g%, leading to an increasing supercurrent responsible for the London field. Finally, we 
remark that in the present case in which only A and A 8 are the nonzero gauge fields, no self-coupling of the London 
field occurs. 
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V. CONCLUSION 



In this paper we have investigated the response of homogeneous superfluid quark matter in thermodynamic equi- 
librium near T c to magnetic fields and rotation. Rotational vortices are topologically stable in the color-flavor locked 
phase, as are magnetic vortices in the isoscalar phase. Table [| summarizes the responses and the corresponding 
corrections to the free energy difference between the superconducting and normal phases. We have ignored several 
aspects, including the modification of the structure of the phases close to T c due to the rotational and magnetic 
energies, effects of the finite strange quark mass m s , and boundary effects. 

The energy correction AG mag due to an external uniform magnetic field H cxt increases the energy of a color-flavor 
locked condensate more than that of an isoscalar condensate, as one can see by c omp aring the color-flavor locked 
result (e 2 /9g|)iJp Xt from Eq. (|60|) and the isoscalar result (e 2 /18gf)H^ xt from Eq. (|l02| ), and noting that o| < 2g\. 
In addition, the energy correction A£2 ro t due to rotation at constant angular velocity u, as given by Eq. (fTq ) in the 
color-flavor locked phase and Eq. (125) in the isoscalar phase, raises the energy of the color-flavor locked state more 
than that of the isoscalar state at temperatures 1 — T/T c > w/gf/i, while being more favorable for the color-flavor 
locked state in a range 1 — T/T c < uj/g^fj,. Although interesting in principle, these terms, for expected magnetic fields 
and rotation rates of neutron stars, affect the phase diagram only in a negligible range around T c and can be ignored. 

The effects of a finite strange quark mass, m s , are also insignificant for the response properties, because a non-zero 
m s does not change the structure of the gradient energy, but merely modifies its magnitude in the color-flavor locked 
phase through the reduction in the Fermi momentum of the strange quarks, and modifies T c . 

Generally, the boundaries of the quark superfluid, e.g., an interface with hadronic matter, contain supercurrents 
induced by the Meissner effect and order-parameter gradients. The resultant energy corrections may possibly compete 
with the bulk corrections listed in Table Q. Magnetic structures have been examined for boundaries of various shapes 
[|lT],Q. On the other hand, the influence of boundaries on rotational structures remains an interesting problem. If a 
rotating neutron star contains hadronic matter with a neutron superfluid and quark matter with a color-flavor locked 
condensate separated by a spherical interface, then the connection of the neutron vortices to the quark vortices is a 
complicated and interesting issue, since the two types of vortices carry a different baryon number and hence a different 
unit of circulation at the interface (3tt//z for quarks and 7t//k for neutrons). For example, a singly quantized vortex 
continuing through the interface would produce a difference in the velocity field around the vortex core, which can 
lead to a Kelvin-Helmholtz instability of the interface 1 33 1 . A connection between three singly quantized neutron 
vortices, each of circulation 7r//x, (or a triply quantized neutron vortex) and a singly quantized quark vortex would 
be free from such velocity differences, since the baryon chemical potential is continuous across the interface. Neutron 
vortices at an interface between superfluid neutron matter and an isoscalar condensate would terminate, while color 
currents responsible for the uniform London field would be arranged in the vicinity of the interface in such a way as 
to prevent the color fields from penetrating into the hadronic region. 
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APPENDIX: HOMOGENEOUS SUPERFLUID 



In this appendix, we summarize the general Ginzburg-Landau analysis for homogeneous superfluid quark matter 
of three massless flavors, temperature T, and baryon chemical potential /x, as developed in I, ignoring the small free 
energy corrections arising from the constraint of overall color neutrality. BCS pairing between two quarks of colors a 
and 6, flavors i and j, spinor indices /i and v, and space-time locations x and y, respectively, is generally characterized 
by a gap matrix, A^ ■ [x — y), which is coupled in the action with the quark spinor [ipt,j {y)\ v and its charge-conjugate 
spinor [ - 0^(x)] M , defined by ip^i = C^ai with C — ij 2 j° (in the Pauli-Dirac representation). The gap A(ai — y) is 
related to the order parameter (tp c (x)'ifj(y)) via the gap equation [Eq. (9) in I]: 
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f j4/_ ..\„ig(x-y) p. (-^ a ) T 



n odd 



x {m c (x)i>(y)})T^\ q , k) + (T[Tp c {x)i> c {y)])T%\q, k)} D$(k - q) , (Al) 

where r^ 11 ) is the full quark-quark-gluon vertex, r^ 21 ) is the full antiquark-quark-gluon vertex, D is the full gluon 
propagator, and the Matsubara frequencies are given by qo = irnvT. (Note that the convention for this order parameter 
is the complex conjugate of the usual condensed matter convention, A ~ (ij)tp).) 

The gap matrix for pairing between quarks with zero total angular momentum, even parity, and aligned chirality 
has the structure [I] 

A(k) = J d{x - y)A(x - y)e lk{x - v) 

= 7 5 [0 + (fco,|k|)A+(k)+0_(fc o ,|k|)A-(k)] , (A2) 
where k is the relative four-momentum of the paired quarks, k = k/|k|, the 

A±(k) = (A3) 

are energy projection operators for noninteracting massless quarks, and ef>± denotes the quark-quark or antiquark- 
antiquark pairing gap. The anti-commutation relations for the quark fields require 

[<f>±(k , |k|)] a6 y = [0±(-fco, \M)]baji ■ (A4) 
In the weak coupling limit, the coefficients a + , /3 + , and /3 + m the Ginzburg-Landau free energy (|l]) are 

a+ = N(^/3)bx(pj , (A5) 



<i = "• = , (A6) 



with the ideal gas density of states at the Fermi surface, 



W3) = i(~) 2 • (A7) 



2tt 2 V3 

Effects beyond weak coupling lead to a no nz ero (3 + through the dependence of the pairing interactions on the gap 



and as well modify all the coefficients ( A5)-(|Aq) mainly through the normal medium corrections to T c and N(fj,/3); 
however, these coefficients remain to be determined in the strong coupling regime. 

For condensates antisymmetric in color and flavor, characterized by Eq. (j33|), substitution of ( |33| ) into Af^°), 
Eq. (0), leads to 

Aft (0) = a\ + {fix + /3 2 T)A 2 (A8) 

with 

a ee 4a+ , /3i ee 16/3+ + 2/3+ , [3 2 = 2/3+ , (A9) 

AEE]T|d a | 2 , T EE^^|d:.d fc | 2 . (A10) 

a ab 

In homogeneous systems, we obtain two types of optimal condensates, as long as the overall fourth-order term in 
Af^ ' is positive. The first, the two-flavor color superconducting (2SC) phase, characterized by T = 1, corresponds 
to an order parameter involving only one component of flavor-antitriplet states, i.e., satisfying 

d R || d G || d B . (All) 
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For T = 1 we focus on the isoscalar, color-antitriplet pairing state with 

( ( l ) +)abij — £abc£ijs(<lc)s ', (A12) 

this state is degenerate in color space under rotation of the color axes. 

The other homogeneous condensation is described by a set of order parameters satisfying T — 1/3 or, equivalently, 

• d G = d* G ■ d B = d* B ■ d R = , \d R \ 2 = \d G \ 2 = \d B \ 2 , (A13) 

corresponding to the color-flavor locked phase Q ; the condensate in this phase is characterized by its symmetry under 
simultaneous exchange of color and flavor. Color-flavor locked states transform into one another under global U(l) 
and flavor (or color) rotation. The simplest among these states is described by (d a )i oc S a i and (d R ) u = (dc)d = 
(dfi) s = kaI the corresponding gap matrix is given by 

{4>+)abij = K A (5 a iSbj - S a jS b i) . (A14) 

For a positive overall fourth order term in the free energy ( |A§| ), the transition from the normal to the superfluid 
state is second order. The resultant state is color-flavor locked for fli > 0, while isoscalar for fa < 0. We remark 
in passing that the critical temperature T c is the same for these two pairing states, which are induced by the same 
instability of the normal phase. Should the coefficient of the fourth order term, /3i + /?2 Y, become negative the phase 
transition becomes first order, as discussed in I. 



[1] B.C. Barrois, Nucl. Phys. B129, 390 (1977). 

[2] D. Bailin and A. Love, Phys. Rep. 107, 325 (1984). 

[3] F. Wilczek, Nucl. Phys. A663, 257c (2000); K. Rajagopal, Nucl. Phys. A661, 150c (1999); K. Rajagopal and F. Wilczek, in 
At the Frontier of Particle Physics, Handbook of QCD, Boris loffe Festschrift, V. 3, edited by M. Shifman (World Scientific, 
Singap ore, 2001), p. 2061 ; D.H. Rischke and R.D. Pisarski, in Proc. Fifth Workshop on QCD, Villefranche-sur-Mer, France 



(2000), [iucl-th/0004016 



[4] M. Alford, K. Rajagopal, and F. Wilczek, Nucl. Phys. B537, 443 (1999). 

[5] N. Evans, J. Hormuzdiar, S.D.H. Hsu, and M. Schwetz, Nucl. Phys. B581, 391 (2000). 

[6] T. Schafer, Nucl. Phys. B575, 269 (2000). 

[7] K. Iida and G. Baym, Phys. Rev. D 63, 074018 (2001), hereafter referred to as I. 

[8] In neutron stars, where the normal flux diffusion time is far greater than the time to cool below the superfluid transition 
temperature ||, the proton superfluid does not have time to expel an external magnetic field, and being a Type II 
superconductor, forms magnetic vortices, even though the magnetic fields are below the lower critical field, H c i , for vortex 
penetration. The expected flux diffusion times here are similarly of order R 2 /c 2 r, where 7? is the neutron star core radius, 
and t is the microscopic scattering time. Thus the flux diffusion times, of order millenia or longer, are again too long 
for magnetic flux expulsion, and one expects that the system must coexist with the fields metastably. In systems with 
color-flavor locked and Type I isoscalar condensates, magnetic fields would be frozen in an intermediate state composed of 
alternating regions of normal material with flux density H c and superconducting material exhibiting incomplete Meissner 
screening, while in a Type II isoscalar condensed system, the fields are frozen in a lattice of vortices. 
[9] G. Baym, C.J. Pethick, and D. Pines, Nature 22 4, 673 (1969). 
[10] D. Blaschke and D. Sedrakian, |nucl-th/0006038 . 

[11] M. Alford, J. Berges, and K. Rajagopal, Nucl. Phys. B571, 269 (2000). 

[12] D.M. Sedrakian, D. Blaschke, K.M. Shahabasyan, and D.N. Voskresensky, Astrofizika 44, 443 (2001). 
[13] M. Alford, J.A. Bowers, and K. Rajagopal, Phys. Rev. D 63, 074016 (2001). 

[14] G. Baym and R.I. Epstein, Astrophys. J. 328, 680 (1988), Astrophys. J. 387, 276 (1992); G. Baym, R.I. Epstein, and B. 

Link, Physica B 178, 1 (1992), Astrophys. J. 403, 285 (1993). 
[15] P.G. de Gennes, Superconductivity of Metals and Alloys, (W.A. Benjamin, New York, 1966). 
[16] B.I. Halperin, T.C. Lubensky, and S.-k. Ma, Phys. Rev. Lett. 32, 292 (1974). 
[17] G. Baym and G. Grinstein, Phys. Rev. D 15, 2897 (1977). 

[18] See, e.g., A.J. Leggett, Rev. Mod. Phys. 47, 331 (1975); P.W. Anderson and W.F. Brinkman, in The Physics of Liquid 

and Solid Helium, Part II, edited by K.H. Bennemann and J.B. Ketterson (Wiley, New York, 1978), p. 177. 
[19] B.D. Josephson, Phys. Lett. 21, 608 (1966). 

[20] K. Iida and G. Baym, Phys. Rev. D 65, 014022 (2002), hereafter referred to as II. 
[21] I. Giannakis and H.-C. Ren, Phys. Rev. D 65, 054017 (2002). 



19 



E.V. Gorbar, Phys. Rev. D 62, 014007 (2000). 

E.J. Yarmchuk, M.J.V. Gordon, and R.E. Packard, Phys. Rev. Lett. 43, 214 (1979). 

R. Blaauwgeers, V.B. Eltsov, M. Krusius, J.J. Ruohio, R. Schanen, and G.E. Volovik, Nature 404, 471 (2000). 

K.W. Madison, F. Chevy, W. Wohlleben, and J. Dalibard, Phys. Rev. Lett. 84, 806 (2000); J.R. Abo-Shaeer, C. Raman, 

J.M. Vogel, and W. Ketterle, Science 292, 476 (2001). 

J. A. Sauls, in Timing Neutron Sta rs, edited by H. O gelman and E.P.J, van den Heuvel (Kluwer, Dordrecht, 1989), p. 457. 



M.M. Forbes and A.R. Zhitnitsky, hep-ph/0109173 



F. London, Superfluids, V. I, (Wiley, New York, 1950). 

G. Baym, in Frontiers and Borderlines in Many-Particle Physics, edited by R. Broglia and J. R. Schrieffer (Soc. Italiana 
di Fisica, Bologna, 1988), p. 330. 

D.H. Rischke, Phys. Rev. D 62, 034007 (2000). 

A. A. Abrikosov, Zh. Eksp. Teor. Fiz. 32, 1442 (1957) [Sov. Phys. JETP 5, 1174 (1957)]. 
W.E. Brown, J.T. Liu, and H.-C. Ren, Phys. Rev. D 61, 114012 (2000); 62, 054016 (2000). 

Lord Kelvin, Mathematical and Physical Papers, V. IV, Hydrodynamics and General Dynamics, (Cambridge, England, 
1910). 



20 



